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Dynamical evolution of heavy quarkonia in a deconfined medium 
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We investigate some consequences of the possibility that heavy quarkonia in a quark-gluon plasma 
possess different (quasi-)bound states, between which transitions are possible. In particular, we show 
that the time-evolution eigenstates in the medium are mixtures of the vacuum eigenstates. This 
leads to abundance ratios of quarkonia that differ from those predicted in statistical models or in 
the sequential-melting picture. 
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I. INTRODUCTION 

Heavy quarkonia provide a convenient testing ground 
for analytical approaches to hadron properties in QCD, 
both in the vacuum and in the presence of a (deconfined) 
medium In the latter case, the original idea of a 
direct link between charmonium suppression and decon- 
finement Q has been refined into more involved predic- 
tions for the behaviors of the various states with rising 
temperature (for recent reviews see Refs. 

Many predictions are formulated in a static picture, 
most noticeably in terms of threshold temperatures, 
above which a given state is entirely "melted", while it 
remains intact below. More dynamical approaches to the 
dissociation and formation or recombination of quarko- 
nia in a medium have been considered in various classical 
kinetic frameworks: a la Boltzmann [^-111: in Langevin or 
Fokker-Planck descriptions [1, [13] or through rate equa- 
tions 

A common feature of these studies is their focus on 
charmonia. This is quite natural, since this corresponds 
to most of the existing experimental results. Now, the 
general expectation is that cc pairs can only exist in a 
quark-gluon plasma (QGP) as either a J or an un- 
bound system. Accordingly, the above mentioned studies 
consider these two possibilities only. On the other hand, 
it is thought that several bottomonia are still bound in 
a QGP just above the deconfinement temperature. This 
opens a richer spectrum of possible behaviors for bb pairs, 
especially if transitions can be induced between different 
bound states. In this Letter, we wish to take this latter 
possibility seriously, and to investigate some of the con- 
sequences of this ansatz, restricting ourselves to inner 
degrees of freedom. 

To this effect, we shall hereafter perform a Gedanken- 
experiment: we put at time i = a bound QQ pair in a 
static, infinitely large QGP at temperature T. Then we 
let the system evolve, and follow the populations of the 
various quarkonium states in time. 

In Section |TT] we introduce our model for the heavy 
quarkonia, the QGP and their interaction, as well as for 
the resulting dynamical equations. Section IIIII contains 
our results for the evolving populations of quarkonium 
states. Further results, in particular involving the exter- 



nal degrees of freedom of quarkonia, will be presented in 
a longer, more technical publication [l^. Eventually, we 
discuss both our model with its underlying assumptions 
and our results in Section HVl 



II. MODEL AND METHOD 

In this Section, we describe our model for the heavy 
quarkonium states and the medium in which they are im- 
mersed, as well as for the coupling between both. Then 
we briefly introduce the equations that govern the dy- 
namics of the QQ-state populations. 



A. Quarkonium in a quark-gluon plasma as a 
{system + reservoir} dissipative system 

Our purpose in the present study is to investigate the 
evolution of heavy quarkonium states under the influ- 
ence of the thermal degrees of freedom of the plasma, 
and especially gluons. Thus, we do not consider (light) 
quarks, which would interact with the heavy quark or 
antiquark through non-thermal gluons: our medium is a 
gluon plasma. 

We model this plasma as a static bath of harmonic 
excitations, whose frequencies span a continuum {wa}- 
The bath is treated as a thermal reservoir, whose ther- 
modynamic properties are not affected by the transitions 
between the various states of the embedded quarkonia. 

The modeling of heavy quarkonia is far more delicate 
than that of the plasma. In this Letter, we wish to iden- 
tify generic behaviors that follow from allowing quark- 
antiquark pairs to be in different (quasi-)bound states in 
a deconfined medium, between which medium-induced 
transitions are permitted. In that view, we deliberately 
adopt an admittedly simplistic quarkonium model, which 
relies on a minimal number of parameters, deferring more 
realistic modeling to further studies. In particular, we 
make a number of assumptions, which we shall further 
discuss in Section HVl 

Our first model assumption is the existence of a non- 
relativistic in-medium quark-antiquark potential (l3 - [l9| 
that admits bound states, thereby neglecting any imagi- 
nary part in the potential [l^. To make the model sim- 
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pie, we consider an attractive Coulomb potential V{r) = 
—a/r, with its well-known spectroscopyo As we want to 
investigate the possible influence of transitions between 
levels, we depart from the exact Coulombic spectroscopy 
for the states that are bound in the vacuum: J/^p, Xc 
and ip' on the one hand, T, Xb, T', x'b ^'^d T" on the 
other one [21] Those states are modeled with Coulomb 
wave functions, yet with the same energy with respect to 
the ground IS'-state as measured in the vacuum. This 
is clearly an approximation, which allows us to lift the 
degeneracy between e.g. 25 and IP states, and thereby 
permit single-gluon transitions between them. 

Obviously, a single gluon cannot induce a transition be- 
tween color-neutral QQ states. We thus further assume 
that our spectrum of color-singlet states is accompanied 
by a parallel spectrum of color-octet states, which for 
the sake of simplicity we shall denote similarly to their 
singlet counterparts. 

Eventually, we have to model the unbound QQ pairs, 
which should normally constitute a continuum. Con- 
sistency with our model for in-medium bound states 
through a real potential implies that the states in that 
continuum have higher energies than the bound states, 
i.e. the latter are always energetically favored, which is 
far from granted. This is even more a problem with the 
Coulomb potential, which admits arbitrarily large bound 
states with high principal quantum number — which in 
the static Debye-screening picture would appear as un- 
bound. To get rid of those states, we fix the dissociation 
threshold by considering bound states of the Coulomb 
potential as describing unbound pairs: 25 and IP states 
(resp. 35* and 2P states) and the higher excited states 
for cc (resp. 66) pairs. A minimal approach to mimic 
the unbound character of such states consists in forbid- 
ding transitions from them back to the bound ones. In 
our computations, we have fully discarded the scattering 
solutions of the Coulomb potential, and considered two 
or three levels of "dissociated" states, with at least two 
states per level, to estimate the error on our result. 

To leave room for the possible "recombination" of 
heavy quark and antiquark into a quarkonium state [2^ , 
we also slightly modified the model by allowing transi- 
tions from the lowest dissociated states to bound ones. 
The plots we present in Section IIIII are for results from 
this variant of our model. Further plots will be shown 
elsewhere jT3| . 

Eventually, we need to specify the coupling between a 
QQ pair and the plasma. We assume dipolar coupling, 
that is, the gluons only interact through their chromo- 
electric field. Incidentally, we need also assume that the 
Bohr frequencies between QQ states are included in the 



^ We take a = 0.4 for cc pairs and a smaller a ~ 0.25 for bb 
pairs, to account for the running of the coupling constant and 
the smaller size of the ground state. 

^ For given S and L quantum numbers, we consider for simplicity 
a single (2L-t-l)-fold degenerated state. 



continuum of bath frequencies, so that transitions be- 
tween states can be induced. 



B. Evolution of the QQ states in the plasma 

Now that we have specified the ingredients of our 
model, we can turn to the time evolution. Further details 
will be given in a longer publication (isj . here we shall 
merely outline the calculation. 

We expect that the state of the quark-antiquark system 
at a given time should be a statistical superposition of 
(vacuum) eigenstates. Then a natural approach is to use 
the master-equation formalism [l^. Chapter 4]. Within 
the decorrelation approximation — i.e., technically, as- 
suming that the density matrix of the whole system fac- 
torizes into the product of the density matrix p'^'^ of the 
QQ pair and that of the plasma at every time — , which 
amounts to considering an expansion up to second order 
in the coupling potential between the quark-antiquark 
pair and the plasma, the populations (diagonal elements) 
of p'^'^ are governed by the coupled Einstein equations 

--T.^.^^f^^{t) + Y.^,^.^p^,^{ti (1) 

where the transition rates Vi^^ between QQ levels fol- 
low from Fermi's golden rule (except for those we set 
to zero, to mimic the continuum, as explained above). 
These rates involve a sum over the states of the QGP, 
weighted by their respective probabilities. This intro- 
duces a dependence of all Vi^k on the plasma tempera- 
ture T. 

In the following Section, we present solutions to these 
evolution equations for the cc and the richer bb systems. 



III. RESULTS 

The Gedankenexperiment we discussed in Section |T] 
amounts to picking out an initial condition at f = for 
Eqs. ([U, for instance pf^^{t — Q) = 1 for the ground \S 
state and for the excited levels, and to solve the coupled 
system. 

A first observation, independent from the initial con- 
dition, follows directly from the structure of the coupled 
evolution equations: the latter are not diagonal. As a 
consequence, the populations corresponding to the vac- 
uum eigenstates of the potential do not constitute an 
eigenstate of the "evolution operator" for the vector of 
populations that can be read off Eqs. ([T]); rather, they 
are linear combinations of the latter. That is, the differ- 
ent quarkonium states do not evolve independently from 
each other, but they are coupled together by the medium. 
For instance, even if T' is constantly either dissociated 
or decaying into Xbi at the same time it is recreated, with 
different rates, through the excitation of Xh or — when 
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FIG. 1. Time dependence of the populations of bound, recombining and unbound states. Left: cc pairs at T = 2Tc\ right: hb 
pairs at T = ST^. Here and in Fig. [3 Te = 170 MeV. 



we allow for that possibility — through recombination of 
an unbounded hh pair. Thus we cannot have the total 
disappearance of a state in the plasma at finite times, as 
implied by static approaches to quarkonium suppression. 

Let r denotes the smaller (in absolute value) of the 
eigenvalues of the evolution operator for the populations. 
The medium-induced mixing of states is such that af- 
ter some transient behavior, which depends on the spe- 
cific choice of initial condition, the populations pf^ of 
all bound states (and of the unbound states that are al- 
lowed to recombine) evolve with the same characteristic 
time scale F"^. We illustrate this in Fig. [TJ which shows 
the time evolution of the populations of cc (left) and bb 
(right) states at respectively 2Tc and ST^. The charac- 
teristic evolution times for these systems at the chosen 
temperatures within our model are respectively 3.4 and 
3.8 fm/c. Quite obviously, for a given system F"-'^ de- 
creases with rising temperature. 

Past the transient regime, the populations of the var- 
ious bound and recombining states are "equilibrated" 
with each other, in the sense that the population ratios 
remain stable. (There is however no strict equilibrium, 
since QQ pairs are consistently lost to non-recombining 
states). These "stationary" population ratios depend on 
the plasma temperature, as shown for hh pairs in Fig. [21 
The stationary ratios differ significantly from those for 
thermally equilibrated levels, which is easily understand- 
able. The thermal ratios are those which make Eqs. ([1]) 
stationary with all transition rates fulfilling the detailed 
balance condition 

F,^fce-^'/^ = Ffc^,e-^'=/^ Vz,fc. (2) 

Here, the condition is obeyed only by bound (T, Xb, T') 
and recombining (x'j,, T") states, but not by the unbound 
ones, from which there are no transition. It is thus nor- 
mal that the resulting stationary ratios diverge from the 
thermal ratios. 



We have checked that when we do not suppress back 
transitions from the unbound states, but set them ac- 
cording to condition ([5]), then the stationary abundance 
ratios are the same as in thermal equilibrium. 

Given the crudeness of our model, which we wish to 
further discuss in the next Section, we have not at- 
tempted at this stage to make predictions for observables 
in real nucleus-nucleus collisions, which would anyway 
necessitate some extra modeling of the plasma kinetics 
as well as accounting for the hadronic phase. 



IV. SUMMARY AND DISCUSSION 

We have examined the dynamics of the populations 
of heavy quarkonium states in a static QGP within the 




FIG. 2. Temperature dependence of the ratios of the popula- 
tions of excited bb states to that of the ground state (T) in the 
stationary regime (see text). Curves: corresponding ratios in 
a thermally equilibrated system. 
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quantum-mechanical master-equation formalism, assum- 
ing that some elements of the quarkonium spectroscopy 
survive in the plasma, in particular that different bound 
states exist, between which transitions can be induced 
by the medium. Under this assumption, we find that the 
bound states are mixed together by the medium, so that 
they all evolve with the same characteristic time scale. 
This differs from the usual picture of sequential melting, 
inasmuch as no bound state can totally disappear while 
others would survive. 

In addition, we find that after some transient regime, 
the population ratios remain stationary, at values de- 
termined by the plasma temperature, yet different from 
the ratios for quarkonium states in thermal equilibrium. 
While the fact that vacuum eigenstates are no longer 
eigenstates in the QGP is model-independent, the ac- 
tual predictions for the abundance ratios depend on the 
model parameters, and could be used to constrain the 
latter from experimental results. 

Let us now discuss our model. Its key ingredient in our 
eyes is the ansatz of medium-induced transitions — with 
a large enough rate compared to the inverse of the plasma 
lifetime — between bound quarkonium states. This is 
the element which couples the evolutions of the various 
bound states together, irrespective of the details of their 
modeling or of the mechanism responsible for the transi- 
tions. 

The latter is important inasmuch as it selects which 
states are coupled together. Here we wanted to con- 
sider, in analogy to quantum optics, transitions induced 
by the gauge bosons. This choice forced us to invoke 
some hazy "spectrum of color-octet states", paralleling 
that of color-singlets — which is far from being granted, 
given the repulsive nature of the color-octet channel. In- 
stead, we could have conjured non-perturbative effects, 
like some kind of soft color interactions [l^] , to instantly 
turn color octets into singlets in the medium. Such ex- 
planations seem to us to be as disputable as our choice 
in the present context. Alternatively, one could think of 
non-color-exchanging processes, as quasi-elastic collisions 
with off-shell quarks or gluons or scattering of photons, 
provided the latter happen with a sufficient rateH 

If one accepts the possibility of medium- induced transi- 
tions between bound states, then the models for medium, 
quarkonia and their interaction that we have used are 
purposely the simplest ones one can think of, yet still 
realistic enough to illustrate some plausible phenomena. 



The orders of magnitude we obtain for the typical time 
scales for the evolution of quarkonium populations are 
reasonable, which justifies our choice a posteriori. The 
model allows us to compute transition rates between 
bound states (which is the reason why we have kept the 
Coulomb wave functions although with "wrong" ener- 
gies), while these are not known in more realistic models, 
since such transitions have not been investigated before. 
For the dissociation widths of the IS' states, one could 
use the known results at leading [1^ or next-to-leading 
order p7j : this would represent an improvement, yet only 
a partial one. 

Note that the dissociation process that we consider, as 
everyone does, is a classical process, in which energy is 
transferred to the QQ pair. This is inherent to the de- 
scription by a real potential similar to the vacuum one. It 
might turn out that a better description of the transition 
from bound QQ state to unbound quark and antiquark in 
a QGP should involve some tunneling through a barrier, 
as studied in hadronic matter in Ref. [28l |. 

The master-equation approach we have used relies on 
a few assumptions, which we shall detail elsewhere }1S] . 
In short, these amount to assuming — as is also done 
in Boltzmann, Langevin or Fokker-Planck formalisms — 
that the typical time scale of plasma correlations is small 
against the characteristic time scale of the QQ-plasma 
interaction, i.e. the formalism implicitly rests on a "weak- 
coupling" assumption. We are investigating alternate ap- 
proaches that do not make use of this hypothesis [2^ , yet 
this seems only feasible at the cost of some alternative 
approximations . 

Eventually, we have assumed dipolar coupling between 
a QQ pair and the plasma, discarding chromomagnetic or 
quadrupolar and higher order chromoelectric couplings. 
This is a large- wavelength approximation — which might 
be disputable for gluons that should resolve the structure 
of the bound states — that can be released when using 
a more realistic model of the quarkonia in the plasma. 
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